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<_ 

QQ ' The method of calculating electroweak decay matrix elements between heavy-heavy and heavy- 

light meson states is developed in the framework of relativistic quark model based on the quasipoten- 
tial approach in quantum field theory. This method is applied for the study of exclusive semileptonic 
^1 B ^ •k{p) decays. It is shown that the large value of the final Tr{p) meson recoil momentum allows 

CN . for the expansion in inverse powers of 6-quark mass of the decay form factors at — 0, where 

is a momentum carried by the lepton pair. This l/rUb expansion considerably simplifies the anal- 
ysis of these decays and is carried out up to the second order. The g^-dependence of the form 
factors is investigated. It is found that the q^-behaviour of the axial form factor Ai is different 
from the other form factors. It is argued that the ratios T{B — > pev)/T{B — > nev) and Tx^/Vt 
are sensitive probes of the Ai g^-dependence, and thus their experimental measurement may dis- 
criminate between different approaches. We find T{B — * -nei/) = (3.0 ±0.6) X \Vub\^ X lO^^s-i and 
V{B — > pev) = (5.4 ± 1.2) x x lO'^^s"'^. The relation between semileptonic and rare radiative 

-B-decays is discussed. 

JlLl' PACS number(s): 12.39Ki, 13.20He 



1 INTRODUCTION 



The investigation of semileptonic decays of B mesons into light mesons is important for the determination 
i of the Cabibbo-Kobayashi-Maskawa matrix element T4b, which is the most poorly studied. At present the 

value of Vub is mainly determined from the endpoint of the lepton spectrum in semileptonic i3-decays Q . 
Unfortunately, the theoretical interpretation of the endpoint region of the lepton spectrum in inclusive 
B Xut-v decays is very complicated and suffers from large uncertainties |^ . The other way to determine 
Vub is to consider exclusive semileptonic decays B —> TT{p)ei'. These are the heavy-to-light transitions 
with a wide kinematic range. In contrast to the heavy-to-heavy transitions, here we can not expand 
matrix elements in the inverse powers of the final quark mass. It is also necessary to mention that the 
final meson has a large recoil momentum almost in the whole kinematical range. Thus the motion of 
final 7r(p) meson should be treated relativistically. If we consider the point of maximum recoil of the 
final meson, we find that 7r(p) bears the large relativistic recoil momentum | Amaa;! of order mb/2 and the 
energy of the same order. Thus at this kinematical point it is possible to expand the matrix element of 
the weak current both in inverse powers of &-quark mass of the initial B meson and in inverse powers of 
the recoil momentum | Amaa,! of the final 7r(p) meson. As a result the expansion in powers l/uib arises for 
the B — > 7r(p) semileptonic form factors at — 0, where is a momentum carried by the lepton pair. 
The aim of this paper is to realize such expansion in the framework of relativistic quark model. We show 
that this expansion considerably simplifies the analysis of exclusive B —^ Tr(p)ev semileptonic decays. 

Our relativistic quark model is based on the quasipotential approach in quantum field theory with 
the specific choice of the qq potential. It provides a consistent scheme for calculation of all relativistic 
corrections at a given order of jc? and allows for the heavy quark l/mq expansion. This model has 
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been applied for the calculations of meson mass spectra Q , radiative decay widths ^ , pseudoscalar decay 
constants heavy-to-heavy semileptonic |^ and nonleptonic Q decay rates. The heavy quark l/mg 
expansion in our model for the heavy-to-heavy semileptonic transitions has been developed in up to 
l/rrig order. The results are in agreement with the model independent predictions of the heavy quark 
effective theory (HQET) [|. The l/mj, expansion of rare radiative decay form factors of B mesons has 
been carried out in jl^ along the same lines as in the present paper. We have briefly presented the results 



for B — > neh' and B —^ pev decays in rcf. 1 11 , where the expansion up to the first order in 1/ mh has 
been carried out. In the present paper we extend the analysis up to the second order and give a detailed 
discussion of the expansion method and results. 

The paper is organized as follows. The relativistic quark model is described in Sect. 2. In Sect. 3 we 
give the detailed description of the method of calculating decay matrix elements between heavy-heavy and 
heavy-light meson states, based on the quasipotential approach. We show that the heavy-to-heavy decay 
matrix elements can be expanded in inverse powers of the heavy quark masses at zero recoil of the final 
meson. On the other hand, the heavy-to-light decay matrix elements can be expanded in inverse powers 
of the initial heavy quark mass at the maximum recoil of the final light meson. These expansions permit 
the calculation of decay matrix elements with the account of relativistic effects. In Sect. 4 the method is 
applied to the calculation of the semileptonic B 'K{p)ev decay form factors. Our numerical results for 
the form factors and decay rates are presented in Sect. 5. Therein we discuss the g^-dependence of the 
form factors and the relations between semileptonic B — > pev and rare radiative B — > p{K*)"f decays. 
Sect. 6 contains our conclusions. The formulae for the form factors at the point of maximum recoil of 
the final light meson are given in Appendix. 



2 RELATIVISTIC QUARK MODEL 

In the quasipotential approach [T^ meson with the mass M and relative momentum of quarks p is 
described by a single-time quasipotential wave function \I'j\/(p), projected onto positive-energy states. 
This wave function satisfies the quasipotential equation 



The quasipotential equation can be transformed into a local Schrodinger-like equation [ p^ 
where the relativistic reduced mass is 



(1) 



(2) 



= 4M5 ' (3) 

and the square of the relative momentum on the mass shell is 

^ ' 4M2 ' ^ ' 

nia.b are the quark masses. While constructing the kernel of this equation V{p, q; M) — the quasipotential 
of quark-antiquark interaction — we have assumed that effective interaction is the sum of the one-gluon 
exchange term with the mixture of long-range vector and scalar linear confining potentials. We have also 
assumed that at large distances quarks acquire universal nonperturbative anomalous chromomagnctic 
moments and thus the vector long-range potential contains the Pauli interaction. The quasipotential is 
defined by §: 

y(p,q,M) = Uaip)ut{-p){^asD^,{khii^^ + KL/(k)r^r,;^ + Vf,„f{k)}ua{q)ut{~q), (5) 
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where as is the QCD coupling constant, D^^, is the gluon propagator; 7^ and u{p) are the Dirac matrices 
and spinors; k = p — q; the effective long-range vector vertex is 

^^(k)=7^ + ^a^,fc^ (6) 

K is the anomalous chromomagnetic quark moment. Vector and scalar confining potentials in the nonrel- 
ativistic limit reduce to 

VZnf{r) = il-e){Ar + B), Vf,,,f{r) ^ e{Ar + B), (7) 

reproducing V^°^l^i{r) = Vf^^j: + V^^j: = Ar+B, where e is the mixing coefficient. The explicit expression 
for the quasipotential with the account of the relativistic corrections of order jc? can be found in ref. 

. All the parameters of our model: quark masses, parameters of linear confining potential A and B, 
mixing coefhcient e and anomalous chromomagnetic quark moment n, were fixed from the analysis of 
meson masses and radiative decays Q. Quark masses: rrih = 4.88 GeV; rric = 1.55 GeV; = 0.50 
GeV; rriu^d = 0.33 GeV and parameters of linear potential: A — 0.18 GeV^; B = —0.30 GeV have 
standard values for quark models. The value of the mixing coefficient of vector and scalar confining 
potentials e = —0.9 has been primarily chosen from the consideration of meson radiative decays, which 
are very sensitive to the Lorentz-structure of the confining potential: the resulting leading relativistic 
corrections coming from vector and scalar potentials have opposite signs for the radiative Ml-decays Q . 
Universal anomalous chromomagnetic moment of quark k — —1 has been fixed from the analysis of the 
fine splitting of heavy quarkonia ^Pj- states j3|. 

Recently we have considered the expansion of the matrix elements of weak heavy quark currents 
between pseudoscalar and vector meson states up to the second order in inverse powers of the heavy 
quark masses [||. It has been found that the general structure of leading, subleading and second order 
1/toq corrections in our relativistic model is in accord with the predictions of HQET. The heavy quark 
symmetry and QCD impose rigid constraints on the parameters of the long-range potential of our model. 
The analysis of the first order corrections H] allowed to fix the value of effective long-range anomalous 
chromomagnetic moment of quarks k = — 1, which coincides with the result, obtained from the mass 
spectra 1^. The mixing parameter of vector and scalar confining potentials has been found from the 
comparison of the second order corrections to be e = — 1. This value is very close to the previous one 
e = —0.9 determined from radiative decays of mesons Therefore, we have got QCD and heavy quark 
symmetry motivation for the choice of the main parameters of our model. The found values of e and 
K imply that confining quark-antiquark potential has predominantly Lorentz-vector structure, while the 
scalar potential is anticonfining and helps to reproduce the initial nonrelativistic potential. 



3 MATRIX ELEMENTS OF ELECTROWEAK CURRENT BE- 
TWEEN HEAVY-HEAVY AND HEAVY-LIGHT MESON 
STATES 

The matrix element of the local current J between bound states in the quasipotential method has the 
form (l|] 

(MV;.(0)|M) = J ^^§^,,(p)r4p,q)*M(q), (8) 

where M(M') is initial (final) meson, r^(p,q) is the two-particle vertex function and are the 

meson wave functions projected onto the positive energy states of quarks. 

This relation is valid for the general structure of the current J = Q'GQ, where G can be an arbitrary 
combination of Dirac matrices. The contributions to T come from Figs. 1 and 2. Thus the vertex functions 
look like 

rW(p,q) =SQ,(pi)GuQ(gi)(27r)3^(p2-q2), (9) 
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and 



L eQ(fci) +eQ(pi) 
+V^(P2 - q2) 7i°G}^Q(gi)u,(g2), (10) 

where ki = pi - A; = qi + A; A = pm-Pa/'; £(p) = (m^ + p2)i/2. 

, ) e{p) ~ (to70 + 7"(7p)) 
^ ^^^^ ^^(^^ ■ 



and 



Pi, 



2 = s^,2{p)^±Y.n('\pM')p' 



i=l 
3 



here 



n 



qi,2 = S,Apf-^±Y^n^^\pM)q\ 

i=l 



:0Mf„)_iM _/ Pi S-. + ^ 



Note that the contribution F^^^ is the consequence of the projection onto the positive-energy states. The 
form of the relativistic corrections resulting from the vertex function F'^' is expHcitly dependent on the 
Lorentz-structure of gg-interaction. 

The general structure of the current matrix element (^) is rather complicated, because it is necessary to 
integrate both with respect to d^p and d^q. The 5-function in the expression for the vertex function 
r(i) permits to perform one of these integrations. As a result the contribution of F^^^ to the current 
matrix element has usual structure and can be calculated without any expansion, if the wave functions 
of initial and final meson are known. The situation with the contribution F^^^ is different. Here instead 
of ^-function we have a complicated structure, containing the potential of gq-interaction in meson. Thus 
in general case we cannot get rid of one of the integrations in the contribution of F^^^ to the matrix 
element (^). Therefore, it is necessary to use some additional considerations. The main idea is to expand 
the vertex function F^^^^ given by ([lO|), in such a way that it will be possible to use the quasipotential 
equation (|l|) in order to perform one of the integrations in the current matrix element ^ . The realization 
of such expansion differs for the cases of heavy-to-heavy and heavy-to-light transitions. 



3.1 Heavy-to-heavy decay matrix elements 

At first we consider the heavy-to-heavy meson decays, such as semileptonic B — > Dev decays, and 
radiative transitions in quarkonia and B* Bj, D* Dj. Here we have two natural expansion 
parameters, which are the heavy quark masses in initial and final meson. The most convenient point for 
the expansion of vertex function F*-^) in inverse powers of the heavy quark masses for semileptonic decays 
is the point of zero recoil of final meson, where A = 0. For radiative decays the momentum transfer 

is fixed |A| = — "2M ■ difference of initial and final meson masses is proportional to the fine or 
hyperfine splitting and thus |A|/Af = o(l/M^), so zero recoil is a good approximation. 

It is easy to see that F^^-* contributes to the current matrix element at first order of 1/toq expansion 
for transitons between mesons consisting from heavy and light quarks {B, D mesons) H] and at second 
order of v/c expansion for mesons consisting from two heavy quarks of the same flavour (quarkonia T, 
J/^*) [Q. We limit our analysis to the consideration of the terms up to the second order in l/mq or v/c 
expansions. We substitute the Dirac matrices G and spinors u in the vertex function F^^-* and consider the 
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cases of Lorentz-scalar and Lorentz-vector (with Pauli term) gg-interaction potential. Then we expand 
r^^) to the desired order and see that it is possible to integrate either with respect to (Pp or (Pq in the 
current matrix element ^ using quasipotential equation (|l|). Performing these integrations and taking 
the sum of the contributions of F^^^ and P^^^ we get the expression for the current matrix element, which 
contains the ordinary mean values between meson wave functions. Thus this matrix element can be easily 
calculated numerically if the meson wave functions are known. The described method has been applied 
to the calculations of hcavy-to-heavy semilcptonic decays in 0, |§] and radiative decays in 0] . 



3.2 Heavy-to-light decay matrix elements 

Now we consider the heavy-to-light meson decays, such as semileptonic B — > TT{p)ei> and rare radiative 
B K*j decays. In these decays the final meson contains only light quarks {u, d, s), thus, in contrast 
to the heavy-to-heavy transitions, we cannot expand matrix elements in inverse powers of the final quark 
mass. The expansion of P*-^-* only in inverse powers of the initial heavy quark mass at A = does not 
allow to use the quasipotential equation for performing one of the integrations in corresponding current 
matrix element (^). However, as it was already mentioned in the introduction, the final light meson has 
the large recoil momentum almost in the whole kinematical range. At the point of maximum recoil of 
final light meson [| the large value of recoil momentum A.max ~ allows for the expansion of decay 

matrix element in l/mg. The contributions to this expansion come both from the inverse powers of 
heavy toq from initial meson and from inverse powers of the recoil momentum | Amaa;| of the final light 
meson. The large value of recoil momentum jAmaa;! permits to neglect in comparison with A^^^ in 
the light quark energy Sq^Q'{p + A) in final meson in the expression for the matrix element originating 
from P''^-'. Such approximation corresponds to omitting terms of the third order in I/uiq expansion 
and is compatible with our analysis, which is carried out up to the second order. It is easy to see that 
we can now perform one of the integrations in the current matrix element (|^) using the quasipotential 
equation as in the case of heavy final meson. As a result we again get the expression for the current 
matrix element, which contains only the ordinary mean values between meson wave functions, but in this 
case at the point of maximum recoil of final light meson. This method has been applied to calculation 
of rare radiative decays of B mesons in ref. pO[ | and in the next section we use it for consideration of 
B — > 7r(p)ej^ semileptonic decays. 



4 B ^ 7T{p)eiy DECAY FORM FACTORS 
4.1 Decay form factors at = 

The form factors of the semileptonic decays B — > Trev and B — s- pev are defined in the standard way as: 
{nip^)\q-/t,b\B{pB)) - f+{q^){PB+P.)^. + f-{q^)iPB-p.)^., (H) 
(Pbp, 6)197^(1 - l')b\B{pB)) = -{Mb + Mp)A,{q^)el + -AM_(e*ps)(ps +Pp)p 



where q = pb — P-k(p)j e is a polarization vector of p meson. In the limit of vanishing lepton mass, the 
form factors /_ and A3 do not contribute to the decay rates and thus will not be considered. 

It is convenient to consider the decay B 7r(p)ei^ in the B meson rest frame. Then the wave function 
of the final 7r(p) meson moving with the recoil momentum A is connected with the wave function at rest 
by the transformation 

*.(P) a(p) = i?y'«A)^y'(i?rA)*.(p)o(p), (13) 



^In the case of rare radiative decays the recoil momentum of final light meson is fixed at the maximum value Amai- 
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where D^/'^[R) is the well-known rotation matrix and i?^ is the Wigner rotation. 

The meson wave functions in the rest frame have been calculated by numerical solution of the quasipo- 
tential equation (|^) |p^ . However, it is more convenient to use analytical expressions for meson wave 
functions. The examination of numerical results for the ground state wave functions of mesons containing 
at least one light quark has shown that they can be well approximated by the Gaussian functions 




with the deviation less than 5%. 
The parameters are 



Pb = 0.41 GeV; = 0.31 GeV. 

Now we apply the method for calculation of decay matrix elements, described in the previous section. 
At the point of maximal recoil of final light meson 

|A_| = -A^. (15) 

we expand the vertex function F'^-* for the Lorentz-scalar and Lorentz-vector (with Pauli term) qq- 
interactions up to the second order in 1/mfc. Then we substitute the vertex functions T^^^ and F^^^ in the 
matrix element (^) and take into account the Lorentz transformation of the final meson wave function 
(p^). Performing one of the integrations in the current matrix element (|^) (using the 5-function in F'^^ 
and the quasipotential equation in the contribution of F^^^ ) we get for the form factors of i? ^ irev and 
B ~> pev decays the following expressions at = point 

/+(0) = /i^^(0) + £/f')(0) + (l-£)/]:('^(0), (16) 

Ai(0) = 4'^(0) + £Af')(0) + (l~e)Aj'('^(0), (17) 

^2(0) = 4^\0) + £Af')(0) + (l-e)ylp'(0), (18) 
V{Q) = l/(i)(0) + eF^(2)(o) + (i_e)v^^(2)(o), (19) 

where f']^\ f^''^^^\ ^1^2' ^12 ''^^ ^ ^'^^ 1/'S''^(2) are given in Appendix, the superscripts "(1)" and 
"(2)" correspond to Figs. 1 and 2, S and V — to the scalar and vector potentials of gg-interaction. 



4.2 l/rrib expansion for decay form factors 

Let us proceed further and for the sake of consistensy carry out the complete expansion of form factors 
(p^-(p^ in inverse powers of &-quark mass. For this expansion we will use some model independent 
results obtained in HQET [g. 

In HQET the mass of B meson has the following expansion in 1 /mi, j9| 

Mb = mj, + A + ^ + O ( ^) , (20) 

where parameter A is the difference between the meson and quark masses in the limit of infinitely heavy 
quark mass. In our model A is equal to the mean value of light quark energy inside the heavy meson 
A — (Sq) B ~ 0.54 GeV Am^ arises from the first-order power corrections to the HQET Lagrangian 
and has the form 

Am| = -Ai-3A2. (21) 
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The parameter Ai results from the mass shift due to the kinetic operator, while A2 parameterizes the 
chromomagnetic interaction ^J. The value of spin-symmetry breaking parameter A2 is related to the 
vector-pseudoscalar mass splitting 

A2 ~ ^(A/|. - Af|) = 0.12 ± 0.01 GeV^. 

The parameter Ai is not directly connected with observable quantities. Theoretical predictions for it vary 
in a wide range: Ai = -0.30 ± 0.30 GeV^ |, |l6|. 

In the limit mg — s- 00, meson wave functions become independent of the flavour of heavy quark. Thus 
the Gaussian parameter (3b in should have the following expansion |^ 

/3b = /3 - ^ + O (^] , (3 w 0.42 GeV, (22) 



nib 



where the second term breaks the flavour symmetry and in our model is equal to w 0.045 GeV^ 



Substituting (20) in (hq) and (fTQ) we get the l/nib expansion of the recoil momentum and the energy 



of final vector meson: 



Now we use the Gaussian approximation for the wave functions (Q). Then shifting the integration 

variable p in (pq)-(f33) by — ^ — \^max, we can factor out the ^max dependence of the meson 

wave function overlap in form factors /+, Ai 2, V ■ The result can be written in the form 



/+(0) = .F+(ALJ exp(-CALj, (24) 

Ai.2(0) - A,2(ALJ exp(-CALj, (25) 

^(0) = V(ALjexp(-CALj, (26) 

where | A^aa;! is given by (pi) and 



.^2 2A^Ai„, / AfB-Af.(p) 



+ /3'(p) ) (E.ip) + Af.(p) )2 Pi ' V Mb + M, 



here 77 = , — and A is equal to the mean value of light quark energy between B and 7r(p) meson 
states: 

A = (sq) w 0.53 GeV, (28) 
Expanding ( p7| ) in powers of 1/ mb we get 



CAi.„^^,(^l-4-;H.j+0^^-,j^ (29) 

We see that the first term in this expansion for the decay into p meson is large. Really, AMp/nib « 0.63. 
The value of this correction is also increased by the exponentiating in (Q)-(^6|). Therefore, we conclude 
that the first order correction in l/mb expansion for the form factors of i? — > T:{p)ev decay, arising from 
the meson wave function overlap, is large.^ Thus, taking into account that our method of calculating decay 
matrix elements does not require the expansion of the meson wave function overlap, we use unexpanded 



expression ( [27D in the exponential of the form factors (p4[)-(26) 

^ The same situation occurs for the rare radiative B decays |lo| . 



7 



In contrast to the meson wave function overlap the factors J^i{Af^^^^), Ai^2(^max) ^i^max)' 
which originate from the vertex functions r^^^'*^^) and Lorentz-transformation ([l3| ) of the final meson 
wave function, have a well defined 1/mf, expansion. The first and second order corrections are small. 
Substituting the Gaussian wave functions (|l^ ) in the expressions for the form factor ( [l6|)- ( |l9|) a nd (^8|)- 
(|69|), with the value of anomalous chromomagnetic quark moment k = — 1, and using "^)^^) and the 
expansions (p0|)-(p3|), we get up to the second order in l/rrii, expansion: 
a.) B ^ TTeiy decay 



max/ 



4^)(ALJ + e^+^'^(A^_) + (1 - e):Fl^'>{A 
n\i + —Xi 

nib 

r^2 



max J 7 



3 \e, 



m„ 



2Ar] - niq - 2M^ 



-A 



.F^(')(ALJ = N\- — {2RZi + {l + 2R)Z2) + —[-=-mgRZ3 + i{p')[l-^R 



(30) 



(31) 



-2mg(A/^-2(£,)) 



1 



3 \ e, 



{A'U{1 -R) + 2m.g{l - 2R) - 2AR) 



1 r 



p^)(l - 2R) + -A77 2(e,)(l - 2R) - AU{1 - R) + 2AR 



l/p2 



3 \ £0 



rrib 



2A-L 



2AU [1-Ir] +4m„ ( 1 - T-R 



-AR 



{1 + 2R) - 2m, (1 - 2R) Z2 



(32) 



(33) 



h) B ^ pev decay 

^i(a: 



•^1 I'^m 



' ) 

max I 



max f 



= 4'^(A^^,J + sAT\Ai,,,) + (1 - e)A^ ^^^(A^,J; 

= iv| 1 + — A:_ + ( r_ + Afp(Mp - raq) + A(Afp + m,) - \Ai]{hmq + 3A/p) 



lS(2) 



y(2), 



(34) 



,2 



3 \ e, 



2Ar/ - - 3A/p 



1 



mi, 

2 

3 \eo 



= AT _(Mp - 2(£,)) + ^ -2(A/p + A + 3m,)(A/p - 2(£,)) 



(Afpi? + mg{l + 3i?) + A(l + R)) + -(p2)(l + 3R) 



A7] A/p(2 + i?) - (e,)(7 + 3i?) + A(3 + R) + 



1 



l/p2 



3 \ e, 



(1 + i?) 



1 



TO6 



iV<^ — (4i?Zi + 2RZ2) + — -{mq - 3Mp)RZ3 - lO(p^) [1--R 



Sn + m„ 



2AL 



8to„ 



+2 2m,(l - i?) - 2A/pi? + r?- 



+2A77 



) A?7 Afp 1 - -i? - A + 1 - —R 
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(35) 



(36) 
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+2(3A - M, - (s,)) -Arj[l-^R\-f/^ 



(37) 



^2 (A 

4^) (A 



max J 



= 4'^(ALJ + sA^'i^'La.) + (1 - £)^r^^(A^,J; 

= iv( 1 + — X_ + ^ ( r_ - Mp{3Mp - m,) + A(Mp - m,) - \kT]{Mp - 7m,) 



triax } 



4('HaLJ = ^r'^(A™a.) + iV^^4Afp(A/p-2(£-,)) 



VfA^ ) 

^ \ max J 

v(i)(a2 ) 



(38) 



nib 



2^2 



- /l5(2) 



3 \eq + niq / \ 



1 



1 /4 



(39) 
(40) 



iV<j — (4i?Zi + 2RZ2) + [ ^{Mp ~ mg)RZ3 - 2(p^) {l-^B. 



£0 + mo 



mi\3 



-Mp (1 - Hi?) + 3 A + 2m,(l + R) 



-2 lmg{l + 'SR)-3MpR + r]-^^R\ Z, 



-2A77 - 



1 \ r 



+3A - Mp - {sq) + Ar/ ( 3 + y i? 



) A77 MJ 1 + -i? + A + m J 3 + — i? 
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1 /p2 



3 \ e. 



= V(i)(ALJ + sV^^^HaLJ + (1 - e)V^(^)(ALJ; 



1 



1 



7 r 



Nl 1 + — X+ + — F- + Afp(A/p - m,) + -TO^ + A(2m, - Afp) + -Kr]M, 



3^^ 



3 \eq + niq 



{Mp - m,) 



(41) 
(42) 

p 

(43) 



V^^'HaLJ = - — (Afp-2(e-,)) + -2 -2(A/p-A-m,)(Afp-2(£-,)) 



nib 



3 \ e„ + mq 



(Afp(l -R) + mq{2 - 3i?) - AR) + -(p^)(2 - 3i?) 



Ar; -Afp(l + 2i?) + 6{£q)R + 2A(1 - R) 



2/p2 



3 \ e. 



(2-i?) 



(44) 



V^(2)(^2^^j = AT — (8i?Zi-2(l-i?)Z2) + — ( --(A^p-m,)i?Z3 + 20(p■ 



mh 



e„ + m„ 



2Afo 
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y 

■2 A 02 



22 - 

-iR] + y A + 4mg (5 + 2R) 



H-2 A+ M,«H-4M 



where A^ 



2fa/3x(p) 



-A^r/^ ( 17 
3/2 



£0 + "^0 



16 



{1- R) + mq{l + R)j Z2 
19 



-Afp + y A + 17m, 



(45) 



/3b ^' 



3/2 



is due to the normaUzation of Gaussian wave functions in 



p2 + m2 + A'^rf, i. e. the energies of Hght quarks in final hght meson acquire additional 
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contribution from the recoil momentum. The averaging (. . .) is taken over the Gaussian wave functions 
of B and 7r(p) mesons, so it can be carried out analytically. For example, 

= / 7^^^(p)(p)£g(p)*s(p) = ^-T^^e^i^i(z), (46) 

where = + A^?7^ and Ki{z) is the modified Bessel function; z = TO^/(2r7/?^(.^^). Analogous 
expressions can be obtained for the other matrix elements in (p0|)-(p5[). 
We have introduced the following notations: 



3 \eq + nig / 2 
1 / P' \ , 1 



3\(e^ + mg)^/ ^'^^ \eg + mq/ 



^3 \eq(£5 +mg)/ \eq + mq) 
Zs = (j^ -2\(A + Af,(p)+3m,) 

\ [Eq + mqY I 



and 



R 



max ) +mb V5 + 2 

5 RESULTS AND DISCUSSION 

Using the parameters of Gaussian wave functions (O) and the value of the mixing coefhcient of vector 



and scalar confining potentials e = — 1 H in the expressions (M), (30)-(|33|) for the B ^ n transition 



form factor /+(0) and the eqs. (|2^), (|2^), (|34|)-(|4^) for the B ^ p transition form factors ^i(O), yl2(0) 
and 1^(0) we get 

f^^'^{0) = 0.20 ±0.02 y^^''(0) = 0.29 ±0.03 

Af^''(0) = 0.26 ±0.03 ^^^''(0) = 0.31 ±0.03. (47) 

The theoretical uncertainty in ( ^ ) results mostly from the approximation of the wave functions by 
Gaussians (^4|) and does not exceed 10% of form factor values. In ref. ijlll we have presented the results 
for B — !■ T:{p)eu decay form factors up to the first order in l/m{, expansion. The found values of form 
factors ||ll[ are very close to (|4^) , this indicates that the second order correctons in 1 /mg are small (less 
than 5% of form factor values) . 



We compare our results (47) for the form factors of i? — > TT{p)ei> decays with the predictions of quark 
models |l^, |l^, QCD sum rules ||l^, ^ and lattice calculations ^ in Table 1. There is an 
agreement between our value of /+~*'^(0) and QCD sum rule and lattice predictions. Our B —> pev form 
factors agree with lattice and QCD sum rule ones |pT| , while they are approximately 1.5 times less than 
QCD sum rule results of refs. ||l^, 

To calculate the B ■n{p) semileptonic decay rates it is necessary to determine the q^-dependence of 
the form factors. Analysing the A^^j^ dependence of the expressions (^-(|69|), (p4[)-(p7|) for the form 



factors /+, Ai, A2 and V ^ wc find that the 5 -dependence of these form factors near q =0 \s given by 
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y(g^) = f; + f; iHv(ALJ, (51) 

^here zz; = ; T+{K.a.) and ^i,2(A^,J, V(A2_) are defined by (|2|)-(|6|), . We 

have introduced the function 



1/2 



which in the hmit of infinitely heavy quarks in the initial and final mesons coinsides with the Isgur-Wise 
function of our model In this limit cqs. (^)-(pT|) reproduce the leading order prediction of HQET 

It is important to note that the form factor Ai in (^9|) has a diff'erent g'^-dependence than the other 
form factors (p8|), (|50|), (pT]). In the quark models it is usually assumed the pole ||l^ or exponential |Q 
q^-behaviour for all form factors. However, the recent QCD sum rule analysis indicates that the form 
factor Ai has q^-dependence different from other form factors |jl^, |2l[. In ||2^ it even decreases with 
the increasing as 

^.(,')-(i-^)a,(0).^^J^(1 + »)A,(0)^ (53) 

Such behaviour corresponds to replacing £,{w) in ( p9[ ) by £,{wmax)- 

We have calculated the decay rates of -B — > TT{p)ev using our form factor values at = and 
the g^-dependence (^)-(^) in the whole kinematical region (model A). We have also used the pole 
dependence for form factors f+{q^), A^iq^), V{q^) and A^{q^) - ^J^^^(l + w)j^^^^ (model B), 

which corresponds to replacing the function ^{w) (p3) by the pole form factor. The results are presented 
in Table 2 in comparison with the quark model 0nill, QCD sum rule | p^ , pO[ and lattice (for B —>■ irev) 
p2| predictions. Lattice accuracy, at present, is not enough to estimate B pev rates |24|. We see that 
our results for the above mentioned models A and B of form factor q^-dependence coincide within errors. 
The ratio of the rates T{B — > pev)/T{B — > irev) is considerably reduced in our model compared to the 
BSW and ISGW [Q models with the simple pole or exponential q^-behaviour of all form factors. 
Meanwhile our prediction for this ratio is in agreement with QCD sum rule results |jl^, The absolute 
values of the rates T{B — > irev) and T{B —> pev) in our model are close to those from QCD sum rules |^ . 
The predictions for the rates with longitudinally and transversely polarized p meson differ considerably 
in these approaches. This is mainly due to different g^-behaviour of Ai (see (^9|), (|5^ ) or pole dominance 
model 1^1). Thus the measurement of the ratios T{B pev)/T{B — > nev) and Tl/Tt should provide 
the test of g^-dependence of Ai and may discriminate between these approaches. 

The differential decay spectra for B — > 7r(p) semileptonic transitions, where x = and Ei is 
the lepton energy, are presented in ref. |ll] (see Fig. 3). 

We can use our results for V and Ai to test the HQET relation ||25|| between the form factors of the 
semileptonic and rare radiative decays of B mesons. Isgur and Wise [ ^5[ have shown that in the limit of 
infinitely heavy fo-quark mass an exact relation connects the form factors V and Ai with the rare radiative 
decay B ^ p^ form factor Fi defined by: 

{p{Pp,e)\m(Tf,^q''PRb\B{pB)) it^uTae*''p^BV''pFi{(i^) 

+ [e;(M| - Af2) _ {e*q){pB +Pp)jG2(g2). (54) 
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This relation is vahd for values sufficiently close to g^„„2. = {Mb — MpY and reads: 

^^('^ ) = 2M-B M^Tm; + ^IM^^'^'' 

It has been argued in |2^, , that in these processes the soft contributions dominate over the hard 
perturbative ones, and thus the Isgur-Wise relations ( p5|) could be extended to the whole range of g^. In 
Ip^ we developed 1 / nib expansion for the rare radiative decay form factor Fi (0) using the same ideas as 
in the present discussion of semileptonic decays. It was shown that Isgur-Wise relation is satisfied 
in our model at leading order of 1 /rrih expansion. The found value of the form factor of rare radiative 
decay -B ^ p7 up to the second order in 1/mi, expansion is ||lO|| 

Ff^''(0) = 0.26 ±0.03. (56) 



Using (p35|) and the values of form factors (47) we find 



Ff^''(0) = 0.27 ±0.03, (57) 

which is in accord with (|5^). Thus wc conclude that 1/mt, and 1/m^ corrections do not break the 
Isgur-Wise relation (|5^) in our model. 

6 CONCLUSIONS 

We have presented in detail the method of the calculation of electroweak decay matrix elements for 
transitions between different meson states, based on the quasipotential approach in quantum field theory. 
It has been shown that the heavy-to-heavy decay matrix element can be expanded in inverse powers of 
the initial and final heavy quark masses at the point of zero recoil of the final meson. On the other hand, 
the heavy-to-light decay matrix element can be expanded in inverse powers of the initial heavy quark 
mass and large recoil momentum of final light meson at the point of maximum recoil of final meson. As 
a result the expansion of the heavy-to-light decay matrix element in inverse powers of initial heavy quark 
mass arises. This method permits the calculation of various radiative and weak decays of heavy mesons 
with the complete account of relativistic effects. 

This method has been applied to the investigation of the semileptonic decays of B mesons into light 
mesons. The recoil momentum of final tt{p) meson is large compared to the 7r(p) mass almost in the 
whole kinematical range. This requires the completely relativistic treatment of these decays. On the 
other hand, the presence of large recoil momentum, which for = is of order TO{,/2, allows for the 
1/mb expansion of weak decay matrix element at this point. The contributions to this expansion come 
both from the heavy 5-quark mass and large recoil momentum of the light final meson. 

We have performed the 1 /rrib expansion of the semileptonic decay form factors at = up to the 
second order. The g^-dependence of the form factors near q^ = has been determined. It has been found 
that the axial form factor Ai has a g^-behaviour different from other form factors (see (|48|)-(|5l|)). This 



is in agreement with recent QCD sum rule results 20, |2l|. The ratios r{B pei>)/T{B — > ireiy) 
and /Tt are very sensitive to the g^-dependence of Ai , and thus their experimental measurement may 
discriminate between different approaches. 

We have considered the relation between semileptonic decay form factors and the rare radiative decay 
form factor , obtained in the limit of the infinitely heavy &-quark. It has been found that in our model 
1/mh corrections do not violate this relation. 
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A APPENDIX: HEAVY-TO-LIGHT SEMILEPTONIC DE- 
CAY FORM FACTORS AT = POINT 

a) -B — > weu decay form factor 



2£g ^ ^ I eg{p + Amax) +mq 



^ / £b(p)+mb r ^ Mb-E^ ^ (pA 



/ 1 1 / 2 2n/^ K-M„ 

+^ ''>[eq{p + Ar,,^,)+mq^2^J)^^^'''^^y\2mb{eq{p + Arnax) + 

1 1 \ Mb - E^ / 1 1 



£q(p)+TOg £q{p+ Ajnax) +mq}^\'^^^^^' ^^^^ 



_|_ (P'^max) 



max 

)+mq £q(A 

max max 

^q'(^max) 



2A2 



-(Mb-£;.) 



^g(^max)(^g('^max) H" 

1 



^max 



max max 

)+mq) mfc(£fc(A 

max 

X (^Mb +M^- Sb{p) - eq{p) - 2eq (^p + 

eq{Amax)-mq Mb - E^ 

2mh{£b{A.max) + mb) 2mb£q{Amax) 
. M^max)-mq f f 2£q 

£q{Amax) + rUg \ ^ \ E^ + M^ 



^ pI+pI 

2{eq{p)+mq) 



E^-M^ 



£q{Amax){£q{p + Amax) + rUq) 



2 



Mb-E^ f 1 1 



E^ + M^ \eq{A 

max max 

) + mq) mb(£b(A 

max ) + nib) 

2£ 

E„ + M^ 

{E, - M,){Mb ~ Skip) - Sq (p) ) Mb-E^ 



X ( Mb + M^ - Ebip) - £q{p) - 2£, p + -= -2— A„,a.T 



jnax 

) + m;,)(£,(A 

max 

) + TTlq) {E^ + M^)mb£q{Amax) 



X 



eq{Amax) - mq f f 2£q \ 

— r— I iW„ Ztq I p -\- — T-^'-^max I 

£q{Amax)+mq\ ^\ E^ + M^ J 



.*s(p), (59) 



K+M, "'"VV 2£,(p + A„„,) 
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pI+pI 



{eq{p) + rUg) \eq{p) + niq eq{Amax) + 



max 

) + ruq V£&(A max )+mb eq{A 

max ) + ruq 



1 



max 

) + rUq eb{A. 

max 

X [Mb +M„- Ebip) - eq{p) - 2e, (p 



E^ + M^ 



+ 



(p'^max) 



1 



{£q{Amax) - rUq) 



£b{Amax) + rrib £q{Amax){£q{Amax) + rUq) 
1 1 



+{Mb - E„) 



+ ■ 



max max ) +mb 

2ea 



x{Mb + M^- Ship) - eq{p) -2eq[p + 
-M 



-A„ 



£q{j)) +mq 



E^ + M^ 

B-M„- Sb{p) - Sqip) + 2£q {p + e'^I^m^ Ama^ 



Mb - 



Mb + M^- Ship) - Sqip) - 2£q (p + ^^^Amax) 



max 

1 



Mb-e, 



^£q{Amax)+mq 2mb \ £q{Amax) + mq 

b) B — > pev decay form factors 

2^MbMp [Sp [ (Pp ^ f 2s, 



Mb " Sbip) - Sqjp) 
£b{A 

max ) + TOfc 



>^b{p), (60) 



A«(0) 



MB + Mp]j Mp 

£b{p) +mb 
2mb 



J (27r)3 ^ 



P + 



Ep + Mp 



l£g{p+A 

max ) + rUq 

2eq(p + Amax) 



1 + 



2mb{£q{p + Amax) + rUq) 



{Ep - Mp) 



Ps + Plj 
£q{p) +mq 



(61) 



Af<')(0) 



2^MbMp Ep f (Pp 



Mb + Mp y Mp J (27r)3 



*P P 



2eo 



£^g(^maa;) ^qi^^max') 

(pl+pl){Ep-Mp) 



Ep + Mp 

Mp - 2eq (p - 



2eo 



Ep + Mp 



l sq{p + Amax) + nig 
2eq{p + Amax) 



-Ar. 



1 



2e„ 



2mb{£q{p) + mq){eq{A 

max 

) + TTlq) [eb{A 

max 

X ^Mb + Mp - Ebip) - Sq{p) - 2eq (p 
1 



(MB-£6(p)-e,(p)) 



+ 



Ep + Mp 

(P^max) £q(Amax) - 



A2 



1 



max ) +mb) 
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Me + Mp- Ship) - Egip) -2eg(p + 



+ 



-{Mb - eb{p) - £g{p)) 



2^MbMp Ep_ f d^p 
Me + MpV Mp 



^*b(p), 
2ea . 



Ep + Mp 



max 
2£g(p + Amax) 



ipl+plKEp-Mp) 



{Sgip) + mg){Sg{Amax) + rUg) \ {E g{p) + TTlg ) (fffe ( A^aj; ) + mfo) 



1 



{eq{A„iax) +mq)' 



(^Mb + Mp- Ship) - eq{p) - 2eq (^p + 



2s„ 



Ep + Mp 



eq{p) + mq 



1 



xlMp~2eq(p + 



1 



2e, . 



Ep + Mp 



(Mb - Ship) - eq{p)) 



1 



max 

) + nib V3(eg(A 

max ) 

X ( - 2£ ( p + „ , A^ 



1 



1 



+ 



2{Mb - Sb{p) - eq{p)) 



/ 2e 

Mb + - Sb{p) - £g{p) -2£g[p+ J L 
1 m„ 



mace max max 

2£o . 



X Ms + Mp - £fe(p) - £g(p) - 2£ J p + 



Ep + Mp 



■*b(p), 



(62) 



(63) 



^(1) _ Mb + Mp2^e;Mp r d^p / 
- 2JM^pEp + MpJ {2^f^'V 



2£„ 



1 + 



Mr 



/£g(p + A 

max 

2£,(p + Ato„^) 

P'. 



Sg{p + Amax) + rUgJ 2mb{eg{p + Amax) + rUg) 



pI+pI 



{eqip + Amax) +mq){eq(j>) +mq) \ 2mb 



Ep + Mp 



Mp 

Mb 



(P^max) ( eq{p + Amax) - mg , ^p f Ep . . , 



(64) 



Me + Mp 2,/e;M'p f d^p 



2/MbMp Ep + Mp J (27r)3 " 



^ . . £g(Am.aa;) "mq _|_ 2 ^"P 



M„\ 1 



2£„ 



Sp + Mp 



mace 

) + mg Mb y £g(A, 



pI+pI 



Mp-2sg(p + 



leq{p + A 

max ) + mq 
2eq{p+ Amax) 

2e„ . 



Ep + Mp 



2mb{sg{p) + mq) 



£6(A 

max ) +mb 



15 



+ 



+ 



{Mb - Ship) - £g{p)) 



+ 



A2 2 



max 
2Sa 



.(A \^rr, I ^MB + Mp-ei,{p)-s,{p)-2e,lp~^ 



A, 



+ — {Mb - £b{p) - £g{p)) }>*S(P), 
Till) 



Ms + Mp 2^EpMp f <Pp 



2^MbMp Ep + Mp J (27r) 

/ pI+pI 



*P P 



^3 



2e„ 



Ep + Mp 



2eq{p+ Ajnax) 



Ep + Mp 



(Sgip) + mg){eq{Amax) + rUq) \ (Sqip) + mq)(Sb{Amax) + m) 



1 



£q{A^f;^x^ + TTlq J \ 
r^2 



[MB + Mp- Ebip) - eq{p) - 2Sq (p + 



2e„ 



Ep + Mp 



£q{p) + niq 



x( MB+Mp- Ebip) - eq{p) -2eq[p + 



2e„ 



Ep + Mp 



1 



-A. 



— \\Mp-2e\p 



max 

) + mb V3(eq(A 

max 

) + TO,) mb 

jP^max) \ 1 f £q{'^max)+mq ^ ITlq 

^liax I £q iP) + TOg V £b{Amax ) + TUb 2£q{Amax ) 

X (^Mb+ Mp - Ebip) - eq{p) - 2eq (p + ^ '^^^ Amax 



2e. 



A 



Ep + Mp 



^q^^max) 



+ 



-[Mp-2eq[p + 



3{Mb - £b{p) - eg{p)) 



2Sa 



Ep + Mp 



>*b(p), 



(65) 



(66) 



Mb + Mp f (fip 



i*p P 



2e„ 



2^MbMp J (27r) V 



Ep + Mp 



2/EpWp 



>£q{p + A 

max ) + rUq 



^ ebip)+mb J ^ ^ Px+Py fsqip+A 

max ) + rUq 

V 2to6 1 
(vA. 

+ 

MB + Mp 



max 

) + TOq y 2eq{p + Amax) 
1 



Ep + Mp\ 2m,b{£q{p) +mq) Sq{p + Amax) + rUq 

1 



jP^max) _ S:q{P + ^max) + mq \ j^^^p) 
^max V 2mb J I 



(67) 



d^p 



2^MbMp J (27r) 
f 1 



3*P(P 



2e„ 



En+Mp 



2y/EpMp £q{p + Amax) + mq 



(^max) 



Mp-2eq[p + 



2e. 



Ep + Mp 



eq{p + A 

max 
^g(^max) 



- 2£g P + 



' A 
2ea 



max 

) + mqV 2eq{p + Amax) 



Px+Pt 

2mbiEp + Mp){sq{p) + rriq) 



Ep + Mp 



-Ar, 
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nib 



max ) 
(P^maa: ) 



/ 2e 

Mb + Mp - e6(p) - £q{p) -2eq[p+ — -^r-pA^aa; 

1 



max max 

2eo 



X ( Mb + Mp - Ebip) - Egip) -2eq[p+ ^ ^"^^ A. 



e,(A 



q X'-^rnax ) 



2m6£g(AmQa;) 



Mp - 2e, p 



2£„ 



Sp + Mp 



>*b(p), 



(68) 



2VMbMp y (2^) 



*P P 



^3 



2£„ 



£;p + Mp 



2^£^pMp l eq{p+ /^,nax) + ™q 



+ A„ 



(iJp + Mp)(£,(p) +TOg) Vlegb) + mg)(eh(A„iaa;) + m;,) 

{MB-eb{p)-eq{p)) 

2£o . 



2mb{eqip) + niq) 
1 



) + mb 

Mb + Mp - £fc(p) - £q(p) - 2£ ( p + 



Ep + Mp 



{P^max) 1 A _ eq{i^„iax) + mq 



^max 

Sqip) +mg \ eb{A 

max ) + rrib 



Mb + Mp - Ebip) - Sqip) -2e[p + 



2£„ 



Ep + Mp 



-A, 



>*s(p), 



(69) 



where the superscripts "(1)" and "(2)" correspond to Figs. 1 and 2, S and V — to the scalar and vector 
potentials of qq-interaction; 



M'B"Kip)_ 



2Mb 

and z-axis is chosen in the direction of A.max- 



Ml + M2, , 
2Mb ' 



(70) 
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TABLE 1 Semileptonic B ^ it and B ^ p decay form factors. 



Ref. 



our results 

& 

[181'^ 



|19| 
1201 
|21| 
|22| 
|23| 



0.20 ±0.02 

0.33 

0.09 
0.26 ±0.02 
0.23 ±0.02 

0.35 ±0.08 



0.26 ±0.03 

0.28 

0.05 
0.5 ±0.1 
0.38 ±0.04 
0.24 ±0.04 
0.24 ±0.12 



0.31 ±0.03 

0.28 

0.02 
0.4 ±0.2 
0.45 ±0.05 

0.27 ±0.80 



0.29 ±0.03 

0.33 

0.27 
0.6 ±0.2 
0.45 ±0.05 
0.28 ±0.06 
0.53 ±0.31 



0.30 ±0.14 ± 0.05 0.22 ±0.05 0.49 ± 0.21 ± 0.05 0.37 ±0.11 



° quark models 
QCD sum rules 
lattice 



TABLE 2 Semileptonic decay rates r{B irev), T{B pev) (x|T4f,P x lO^^s ^) and the ratio of the 
rates for longitudinally and transversely polarized p meson. 



Ref. 



T{B^Tiev) T{B->pev) Tl/Tt 



our results 
model A 
model B 

[18|" 
|19|^ 
|20|^ 
1221'= 



3.0 ±0.6 
2.9 ±0.6 

7.4 
2.1 

5.1 ± 1.1 
3.6 ±0.6 

8±4 



5.4 ± 1.2 

5.0 ± 1.2 
26 
8.3 

12 ±4 

5.1 ± 1.0 



0.5 ±0.3 
0.5 ±0.3 
1.34 
0.75 
0.06 ±0.02 
0.13 ±0.08 



" quark models 
^ QCD sum rules 
lattice 



FIGURE CAPTIONS 

FIGURE 1 Lowest order vertex function 

FIGURE 2 Vertex function with the account of the quark interaction. Dashed line corresponds to the 
effective potential (j^). Bold line denotes the negative-energy part of the quark propagator. 
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